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Simulation  of  Reaction  Dynamics  for  Synthesis  of 
Energetic  Materials  and  Resistant  Coatings 


I.  Introduction 

The  objective  of  the  research  supported  by  this  grant  was  the  development  of 
theoretical  and  computational  methods  to  guide  the  efficient  synthesis  of  HEDM  (High 
Energy  Density  Matter)  for  use  in  advanced  propellants  and  POSS  (polyhedral 
oligomeric  silsesquioxanes)  for  use  in  coatings  resistant  to  extreme  conditions  such  as 
heat  and  abrasion.  The  research  centered  on  the  design  of  new  methodology  for  the 
simulation  of  hydrogen  transfer  reactions.  The  projects  included  the  development  of  the 
following  three  types  of  approaches:  grid  methods  for  the  calculation  of  hydrogen 
vibrational  wavefunctions,  nonadiabatic  molecular  dynamics  methods  for  the  simulation 
of  proton  transfer  in  liquids,  and  the  nuclear-electronic  orbital  (NEO)  method  for  the 
incorporation  of  nuclear  quantum  effects  in  electronic  structure  calculations. 

II.  Methods  and  Results 

A.  Grid  methods  for  hydrogen  vibrational  wavefunctions 

The  goal  of  this  project  was  to  develop  grid-based  methodology  for  the 
calculation  of  hydrogen  vibrational  wavefunctions  for  systems  in  which  only  one  or  a  few 
specified  hydrogen  nuclei  are  treated  quantum  mechanically.'’^  We  have  combined  these 
grid-based  methods  with  electronic  structure  calculations^  and  with  mixed 
quantum/classical  molecular  dynamics  simulations.^  ^  These  approaches  are  applicable  to 
the  hydrogen  transfer  reactions  relevant  to  the  synthesis  of  HEDM  and  POSS. 

The  grid-based  methodology  we  developed  to  calculate  hydrogen  vibrational 
wavefunctions  is  called  the  Fourier  Grid  Hamiltonian  Multiconfigurational  Self- 
Consistent-Field  (FGH-MCSCF)  method.’  In  FGH-MCSCF,  the  potential  energy  surface 
for  the  transferring  hydrogen(s)  is  calculated  on  a  grid  (i.e.,  all  nuclei  except  the 
transferring  hydrogen(s)  are  fixed,  and  the  energy  is  calculated  for  the  hydrogen(s) 
positioned  at  each  point  on  the  grid).  The  hydrogen  vibrational  wavefunctions  are 
calculated  by  numerically  solving  the  time-independent  Schrddinger  equation  for  the 
hydrogen  nucleus  (or  nuclei)  moving  on  this  potential  energy  surface.  We  implement  an 
MCSCF  approach  for  the  calculation  of  multidimensional  wavefunctions.  In  this 
MCSCF  approach,  the  hydrogen  vibrational  wavefunction  is  expressed  as  a  linear 
combination  of  single  configurations  that  are  products  of  one-dimensional  wavefunctions 
represented  directly  on  the  grid.  The  variational  method  is  utilized  to  minimize  the  total 
energy  with  respect  to  both  the  configuration  coefficients  and  the  one-dimensional 
wavefunctions.  Since  only  a  relatively  small  number  of  configurations  is  required  for  the 
calculation  of  accurate  hydrogen  vibrational  wavefunctions,  this  MCSCF  approach 
avoids  the  diagonalization  of  large  matrices  required  for  a  full  configuration  interaction 
approach. 

Although  the  FGH-MCSCF  method  decreases  the  computational  expense  of 
calculating  the  vibrational  wavefunctions,  this  grid  method  requires  the  calculation  of  the 
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potential  energy  and  forces  at  each  grid  point.  The  generation  of  the  grid  potential  for 
each  molecular  dynamics  time  step  is  often  the  bottleneck  of  mixed  quantum-classical 
molecular  dynamics  simulations.  To  decrease  the  computational  expense  of  this  step,  we 
have  developed  a  partial  multidimensional  grid  generation  method.^  This  method 
substantially  decreases  the  number  of  potential  energy  calculations  (typically  by  more 
than  an  order  of  magnitude)  by  avoiding  these  calculations  for  grid  points  with  high 
potential  energy.  The  combination  of  the  FGH-MCSCF  and  partial  multidimensional 
grid  generation  methods  has  been  found  to  accurately  describe  ground  and  excited  state 
hydrogen  vibrational  wavefunctions  in  a  computationally  practical  manner. 

We  have  used  this  methodology  in  two  different  contexts.  First,  we  have 
combined  electronic  structure  and  hydrogen  vibrational  wavefunction  methods  to 
calculate  the  hydrogen  potential  energy  surfaces  and  vibrational  wavefunctions  for 
structures  along  a  reaction  path  obtained  from  electronic  structure  methods.  These 
calculations  provide  useful  information  about  the  fundamental  nature  of  the  nuclear 
quantum  effects  for  hydrogen  transfer  reactions.  In  addition,  we  have  utilized  these  grid- 
based  methods  for  mixed  quantum/classical  molecular  dynamics  simulations  of  hydrogen 
transfer  reactions  in  condensed  phase  systems.  Specifically,  we  have  studied  hydride 
transfer  in  liver  alcohol  dehydrogenase  and  dihydrofolate  reductase  enzymes. 

B.  Nonadiabatic  molecular  dynamics  methods  for  proton  transfer  in  solution 

The  goal  of  this  project  was  to  develop  nonadiabatic  mixed  quantum/classical 
molecular  dynamics  methods  for  the  simulation  of  proton  transfer  in  solution.  The 
emphasis  was  on  the  treatment  of  both  the  hydrogen  and  the  donor-acceptor  vibrational 
motion  quantum  mechanically  within  the  framework  of  the  molecular  dynamics  with 
quantum  transitions  (MDQT)  surface  hopping  approach.*’^  We  applied  this  method  to  the 
model  proton  transfer  reaction  AH  —  B  ^  A  —  *HB  in  liquid  methyl  chloride.  This 
project  provides  the  foundation  for  the  simulation  of  hydrogen  transfer  reactions  required 
for  the  synthesis  of  HEDM  and  POSS. 

In  mixed  quantum/classical  molecular  dynamics  methods,  one  or  a  few  nuclei  are 
treated  quantum  mechanically  while  the  remaining  nuclei  are  treated  classically.  The 
MDQT  surface  hopping  method®’^  incorporates  nonadiabatic  transitions  among  the 
adiabatic  vibrational  states.  The  fundamental  principle  of  MDQT  is  that  an  ensemble  of 
trajectories  is  propagated,  and  each  trajectory  moves  classically  on  a  single  adiabatic 
surface  except  for  instantaneous  transitions  among  the  adiabatic  states.  The  adiabatic 
vibrational  states  are  calculated  at  each  classical  molecular  dynamics  step  by  solving  the 
time-independent  Schrodinger  equation  for  the  quantum  nuclei.  The  classical  nuclei 
evolve  on  the  occupied  vibrational  surface  according  to  Newton’s  classical  equations  of 
motion.  Nonadiabatic  transitions  among  these  surfaces  are  incorporated  according  to  a 
probabilistic  algorithm  that  correctly  apportions  trajectories  among  the  adiabatic  states 
according  to  the  quantum  probabilities  (neglecting  difficulties  with  classically  forbidden 
transitions).  These  quantum  probabilities  are  determined  by  integration  of  the  time- 
dependent  Schrddinger  equation  simultaneously  with  the  classical  equations  of  motion. 
The  MDQT  method  has  been  applied  to  proton  transfer  in  solution®  and  in  enzymes.  ‘  In 
these  previous  applications  of  MDQT,  only  the  transferring  hydrogen  motion  was  treated 
quantum  mechanically. 


3 


Sharon  Hammes-Schiffer 
Final  Technical  Report 


AFOSR  Grant  No.  F49620-0 1-1-0046 
Dec  2000-Nov  2003 


During  the  past  grant  period,  we  have  developed  the  methodology  for  treating 
both  the  hydrogen  motion  and  the  donor-acceptor  vibrational  motion  quantum 
mechanically  within  the  framework  of  MDQT.’°  We  have  applied  this  methodology  to 
proton  transfer  in  a  linear  AHB  complex.  In  this  case,  the  treatment  of  only  the  hydrogen 
motion  quantum  mechanically  requires  the  calculation  of  one-dimensional 
wavefunctions,  whereas  the  treatment  of  both  the  hydrogen  motion  and  the  donor- 
acceptor  vibrational  motion  quantum  mechanically  requires  the  calculation  of  two- 
dimensional  wavefunctions.  We  denote  the  former  treatment  as  a  ID  and  the  latter 
treatment  as  a  2D  mixed  quantum/classical  approach. 

In  the  2D  mixed  quantum/classical  method,  the  vibrational  wavefunctions  depend 
on  the  proton  coordinate  r  and  the  AB  vibrational  coordinate  R.  For  each  classical 
configuration  obtained  during  the  molecular  dynamics  simulation,  the  two-dimensional 
time-independent  Schrodinger  equation  is  solved  numerically  on  a  two-dimensional  grid 
using  the  grid-based  methods  described  above.  The  i?-grid  is  generated  by  moving  A  and 
B  along  the  AB  axis  while  maintaining  the  same  orientation  and  center  of  mass.  The  r- 
grid  is  generated  by  moving  H  along  the  AB  axis.  The  translational  and  rotational 
degrees  of  freedom  of  the  AHB  complex,  as  well  as  all  solvent  degrees  of  freedom,  are 
treated  classically.  For  2D  MDQT,  the  nonadiabatic  transitions  are  incorporated  among 
the  two-dimensional  vibrational  states. 

We  have  applied  the  2D  MDQT  approach  to  the  intramolecular 
reaction  AH-B^A~-  in  liquid  methyl  chloride,  where  AHB  represents  a  phenol- 
amine  complex. This  model  was  initially  constructed  by  Azzouz  and  Borgis."  The  aim 
of  our  study  was  to  compare  the  results  from  ID  and  2D  MDQT  to  elucidate  the 
fundamental  issues  arising  from  a  quantum  mechanical  treatment  of  the  donor-acceptor 
vibrational  motion  as  well  as  the  hydrogen  motion.  We  observed  significant  differences 
in  the  behavior  of  the  systems,  particularly  in  the  character  of  the  excited  state  vibrational 
wavefunctions.  Despite  these  differences,  however,  the  calculated  rates  and  kinetic 
isotope  effects  were  qualitatively  similar  for  ID  and  2D  MDQT.  Thus,  these  calculations 
indicate  that  a  classical  treatment  of  the  donor-acceptor  mode  is  reasonable  for  these 
types  of  reactions. 

C.  Nuclear-electronic  orbital  (NEO)  method 

The  goal  of  this  project  was  to  develop  the  NEO  method  for  the  incorporation  of 
nuclear  quantum  effects  into  electronic  molecular  orbital  calculations.’^  In  this  method, 
specified  nuclei  are  treated  quantum  mechanically  on  the  same  level  as  the  electrons. 
Both  electronic  and  nuclear  molecular  orbitals  are  expressed  as  linear  combinations  of 
Gaussian  basis  functions,  and  the  variational  method  is  utilized  to  minimize  the  energy 
with  respect  to  all  molecular  orbitals,  as  well  as  the  centers  of  the  nuclear  basis  functions. 
Significant  correlation  effects  are  included  using  a  multiconfigurational  self-consistent- 
field  method.  We  have  applied  the  NEO  approach  to  the  proton  transfer  reactions 
required  for  the  synthesis  of  POSS.  This  application  has  illustrated  that  inclusion  of 
nuclear  quantum  effects  significantly  influences  the  energetics  of  these  hydrogen  transfer 
reactions. 

The  NEO  approach  may  be  used  to  calculate  structures,  energies,  minimum 
energy  paths,  and  direct  dynamics  for  chemical  reactions.  Typically  only  a  subset  of  the 
nuclei  is  treated  quantum  mechanically.  For  example,  often  only  the  hydrogen  nuclei  are 
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chosen  for  the  quantum  subsystem,  while  all  other  nuclei  are  included  in  the  classical 
subsystem.  The  specified  quantum  nuclei  are  treated  analogously  to  the  electrons,  and 
properties  of  the  system  are  averaged  over  the  mixed  nuclear-electronic  wavefunction. 
Structures  are  optimized  with  respect  to  only  the  classical  nuclei,  and  the  reaction 
coordinate  along  the  minimum  energy  path  is  comprised  of  only  classical  nuclei. 
Similarly,  for  dynamical  calculations  the  potential  energy  surface  depends  explicitly  on 
only  the  classical  nuclei. 

The  NEO  approach  is  particularly  useful  for  the  description  of  hydrogen  transfer 
reactions.  In  this  case,  the  transferring  hydrogen  nucleus  is  treated  quantum 
mechanically,  and  the  reaction  coordinate  of  the  minimum  energy  path  reflects  changes 
in  the  environment  created  by  the  classical  nuclei.  This  is  analogous  to  the  collective 
solvent  coordinate  in  Marcus  theory  for  electron  transfer  reactions.  The  minimum  energy 
path  includes  nuclear  quantum  effects  such  as  zero  point  energy  and  hydrogen  tunneling. 
This  approach  maintains  the  conceptual  picture  of  a  reaction  path,  while  eliminating  the 
difficulties  associated  with  large  curvature  of  the  reaction  path  for  the  transfer  of  light 
nuclei. 

The  NEO  approach  possesses  many  advantageous  attributes.  One  important 
strength  of  the  NEO  approach  is  that  the  nuclear  quantum  effects  are  incorporated  during 
the  electronic  structure  calculation,  rather  than  subsequently  calculated  as  a  correction 
factor.  In  addition,  the  NEO  approach  does  not  invoke  the  Bom-Oppenheimer  separation 
of  the  electrons  and  the  quantum  nuclei.  In  the  multiconfigurational  NEO  approach, 
excited  vibrational-electronic  states  may  be  calculated,  and  nonadiabatic  effects  may  be 
included  in  dynamical  calculations.  Finally,  the  NEO  approach  is  computationally 
practical  for  a  wide  range  of  chemical  reactions,  and  its  accuracy  may  be  improved 
systematically  through  increased  basis  set  size  and  inclusion  of  additional  correlation 
with  larger  numbers  of  configurations  and  extended  methodology. 

We  are  designing  the  NEO  approach  to  study  chemical  reactions.  As  we  have 
shown,'^'''^  small  basis  sets  and  single  configurational  wavefunctions  are  not  adequate  for 
describing  many  chemical  reactions,  particularly  those  involving  hydrogen  transfer,  and 
are  limited  to  the  calculation  of  only  ground  vibronic  states.  In  the  NEO  approach,  the 
nuclear  basis  set  is  comprised  of  ^-type,  p-type,  and  d-type  Gaussian  basis  functions,  and 
multiple  basis  function  centers  for  each  quantum  hydrogen  nucleus  are  allowed. 
Correlation  among  electrons  and  nuclei  is  included  with  a  procedure  analogous  to  the 
complete  active  space  self-consistent-field  (CASSCF)  electronic  structure  method.  The 
combination  of  multiple  basis  function  centers  and  multiconfigurational  wavefunctions  is 
required  to  accurately  describe  the  delocalized,  bilobal  character  of  the  vibrational 
wavefunctions  involved  in  hydrogen  transfer. We  have  illustrated  this  requirement 
through  an  application  of  NEO  to  malonaldehyde.'^’’'*  We  have  also  benchmarked  the 
NEO  method  through  applications  to  H2  and  HF  and  have  obtained  excellent  agreement 
of  the  vibrational  energy  splittings  with  the  experimental  values. 

In  addition,  we  have  developed  the  methodology  to  perform  a  NEO  vibrational 
analysis  and  have  explored  the  physical  implications  of  such  an  analysis. This  analysis 
is  required  for  the  investigation  of  chemical  reactions  on  the  vibrational-electronic 
potential  energy  surface.  In  the  NEO  approach,  the  centers  of  the  nuclear  basis  functions 
are  optimized  variationally  on  the  same  level  as  the  coefficients  in  the  molecular  orbitals, 
and  the  NEO  potential  energy  surface  depends  explicitly  on  only  the  coordinates  of  the 
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classical  nuclei.  The  NEO  vibrational  analysis  involves  the  calculation,  projection,  and 
diagonalization  of  a  numerical  Hessian  to  determine  the  harmonic  vibrational  frequencies 
corresponding  to  the  classical  nuclei.  The  analysis  of  these  frequencies  allows  the 
characterization  of  stationary  points  on  the  NEO  potential  energy  surface  (i.e.,  distinction 
among  minima,  transition  states,  and  higher-order  saddle  points).  This  analysis  also 
enables  the  calculation  of  zero  point  energy  corrections  and  thermodynamic  quantities 
such  as  enthalpy,  entropy,  and  free  energy  changes  for  reactions  on  the  NEO  potential 
energy  surface.  Furthermore,  the  numerical  Hessian  may  be  used  in  the  generation  of 
minimum  energy  paths,  or  intrinsic  reaction  coordinates,  on  the  NEO  potential  energy 
surface.  We  have  applied  this  vibrational  analysis  to  a  series  of  representative  molecules, 
including  HCN,  the  protonated  water  dimer,  and  triazene,  as  well  as  to  the  reaction 
profile  for  an  identity  nucleophilic  substitution  reaction. 

The  NEO  approach  has  been  implemented  in  the  GAMESS  electronic  structure 
program.'^  A  summary  of  the  current  capabilities  of  our  NEO  code  is  as  follows: 

1.  Utilize  the  DZSPDN  nuclear  basis  set  with  the  option  of  multiple  basis  function 
centers. 

2.  Calculate  expectation  values  of  hydrogen  coordinates  with  respect  to  molecular 
orbitals. 

3.  General  options:  conventional  or  direct  calculation  of  integrals,  RHF  or  ROHF 
treatment  of  electrons,  and  the  possibility  of  calculations  with  symmetry. 

4.  Treat  the  quantum  nuclei  as  fermions,  bosons,  or  a  mixture  of  fermions  and 
bosons. 

5.  Calculate  energies  and  gradients  at  the  NEO-HF  level. 

6.  Calculate,  project,  and  diagonalize  the  numerical  Hessian. 

7.  Locate  and  characterize  geometry  stationary  points  (i.e.,  minima  and  transition 
states). 

8.  Calculate  zero  point  energy  corrections  and  thermodynamic  properties. 

9.  Calculate  energies  at  the  NEO-CI  level. 

10.  Calculate  energies  at  the  NEO-MCSCF  level  and  with  state-averaged  NEO- 
MCSCF. 

1 1 .  Calculate  energies  at  the  NEO-MP2  level. 

III.  Summary 

Thus,  we  have  developed  three  different  types  of  approaches  for  the  simulation  of 
hydrogen  transfer  reactions.  The  grid-based  methodology  enables  the  efficient 
calculation  of  hydrogen  vibrational  wavefunctions  in  conjunction  with  electronic 
structure  calculations  and  mixed  quantum/classical  molecular  dynamics  simulations.  The 
nonadiabatic  mixed  quantum/classical  molecular  dynamics  method  allows  the  quantum 
mechanical  treatment  of  both  the  hydrogen  and  the  donor-acceptor  vibrational  motion  for 
the  simulation  of  proton  transfer  in  solution.  The  nuclear-electronic  orbital  (NEO) 
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method  incorporates  nuclear  quantum  effects  in  electronic  molecular  orbital  calculations 
and  enables  the  calculation  of  structures,  energies,  minimum  energy  paths,  and  direct 
dynamics  for  chemical  reactions. 

The  application  of  these  approaches  to  hydrogen  transfer  reactions  required  for 
the  synthesis  of  POSS  and  HEDM  is  providing  information  that  should  aid  in  the  efficient 
synthesis  of  these  materials.  These  calculations  elucidate  the  role  of  nuclear  quantum 
effects  such  as  zero  point  energy  and  hydrogen  tunneling.  Moreover,  they  provide 
insight  into  the  detailed  mechanisms,  as  well  as  the  influence  of  solvent  and  substituents 
on  the  rates  and  yields.  Such  mechanistic  information  is  critical  for  the  efficient 
synthesis  of  POSS  and  HEDM  and  the  production  of  new  functional  species.  POSS  and 
HEDM  are  of  interest  to  the  Air  Force  due  to  the  wide  range  of  technological 
applications,  including  advanced  propellants  and  coatings  resistant  to  extreme  conditions 
such  as  heat  and  abrasion. 
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“Hybrid  Quantum-Classical  Molecular  Dynamics  Studies  of  the  Relation  between 
Enzyme  Motion  and  Activity” 
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Activity” 

19.  Gordon  Research  Conference  on  Computational  Chemistry,  New  London,  New 
Hampshire,  June  30- July  5,  2002  (invited  talk):  “Hybrid  Quantum-Classical 
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(invited  talk):  “Multiconfigurational  Nuclear-Electronic  Orbital  (NEO)  Approach: 
Including  Nuclear  Quantum  Effects  in  Electronic  Structure  Calculations” 
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Biological  Enzymes  Work?” 
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August  18-22,  2002  (invited  talk):  “Theoretical  Studies  of  Proton-Coupled 
Electron  Transfer  Reactions  in  Solution” 
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26.  CIMMS/CALTECH  Workshop  entitled  Molecular  Modeling  and  Computation: 
Perspectives  and  Challenges,  Caltech,  Pasadena,  California,  November  15-16, 
2002  (invited  talk):  “Hybrid  Quantum-Classical  Molecular  Dynamics  of 
Hydrogen  Transfer  Reactions  in  Enzymes” 

27.  Columbia  University,  New  York,  New  York,  January  31,  2003  (biophysics 
seminar):  “The  Impact  of  Enzyme  Motion  on  Activity” 

28.  Sanibel  Symposium,  St.  Augustine,  Florida,  February  22-March  1,  2003  (invited 
talk):  “Hybrid  Quantum-Classical  Molecular  Dynamics  of  Hydrogen  Transfer 
Reactions  in  Enzymes” 

29.  American  Chemical  Society  National  Meeting,  Symposium  on  New  Electronic 
Structure  Methods:  From  Molecules  to  Materials,  New  Orleans,  Louisiana,  March 
22-26,  2003  (invited  talk):  “Incorporation  of  Nuclear  Quantum  Effects  in 
Electronic  Structure  Calculations:  Multiconfigurational  Nuclear-Electronic 
Orbital  Method” 

30.  American  Chemical  Society  National  Meeting,  Symposium  on  Integrating 
Diverse  Computational  Approaches  to  Complex  Problem  Solving,  New  Orleans, 
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Louisiana,  March  22-26,  2003  (invited  talk):  “Hybrid  Quantum-Classical 
Molecular  Dynamics  of  Hydrogen  Transfer  Reactions  in  Enzymes” 

31.  American  Chemical  Society  National  Meeting,  Symposium  on  The  Cutting  Edge: 
Use  of  Computers  in  Teaching  and  Learning  Chemistry,  New  Orleans,  Louisiana, 
March  22-26,  2003  (invited  talk):  “Utilization  of  Computer  Movies  to  Illustrate 
Quantum  Effects  and  Motion  in  Enzyme  Reactions” 

32.  American  Society  of  Biochemistry  and  Molecular  Biology  National  Meeting, 
Session  on  Fundamental  and  Emerging  Issues  in  Enzymatic  Catalysis,  San  Diego, 
California,  April  12-16,  2003  (invited  talk):  “Impact  of  Enzyme  Motion  on 
Activity” 

33.  AFOSR  Molecular  Dynamics  Contractor’s  Meeting,  San  Diego,  California,  May 
18-20,  2003  (invited  talk):  “Nuclear  Quantum  Effects  in  Hydrogen  Transfer 
Reactions:  Polyhedral  Oligomeric  Silsesquioxanes  and  Ionic  Liquids” 

34.  Gordon  Research  Conference  on  Photosynthesis,  New  Hampshire,  June  22-26, 
2003  (invited  talk):  “Coupling  of  Electrons  and  Protons  to  the  Environment” 

35.  Workshop  entitled  Radicals  in  the  Rockies,  Telluride,  Colorado,  July  6-12,  2003 
(invited  talk):  “Proton-Coupled  Electron  Transfer  in  Solution  and  Enzymes” 

36.  International  meeting  entitled  Multidimensional  Quantum  Reaction  Dynamics, 
Freie  Universitat,  Berlin,  Germany,  July  16-18,  2003  (invited  talk):  “Hybrid 
Quantum-Classical  Calculations  of  Hydrogen  Transfer  Reactions” 

37.  Symposium  entitled  Computational  Modelling  of  Catalysis,  Max  Planck  Institute, 
Muelheim,  Germany,  July  16-18,  2003  (invited  talk):  “Hybrid  Quantum-Classical 
Molecular  Dynamics  of  Hydrogen  Transfer  Reactions  in  Enzymes” 

38.  Conference  entitled  Excited  State  Processes  in  Electronic  and  Bio  Nanomaterials, 
Los  Alamos,  New  Mexico,  August  11-16,  2003  (invited  talk):  “Proton-Coupled 
Electron  Transfer  Reactions” 

39.  American  Chemical  Society  National  Meeting,  Symposium  on  Making  and 
Breaking  Chemical  Bonds  in  Gas  and  Condensed  Phases:  Theory  and 
Applications,  New  York,  New  York,  September  7-11,  2003  (invited  talk): 
“Investigation  of  Hydrogen  Transfer  Reactions  with  the  Multiconfigurational 
Nuclear-Electronic  Orbital  Method” 

40.  Central  Regional  American  Chemical  Society  Meeting,  Pittsburgh,  Pennsylvania, 
October  19-23,  2003  (invited  talk):  “Hybrid  Quantum-Classical  Molecular 
Dynamics  of  Hydrogen  Transfer  Reactions  in  Enzymes” 

41.  University  of  Iowa,  Iowa  City,  Iowa,  November  6,  2003  (colloquium):  “Impact  of 
Enzyme  Motion  on  Activity” 
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